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~£T)~. The scattering of two charged strongly interacting particles is described on the basis 

§ ' of the P-matrix approach. In the P matrix, it is proposed to isolate explicitly the 

(N '■ 

background term corresponding to purely Coulomb interaction, whereby it becomes 
possible to improve convergence of the expansions used and to obtain a correct asymp- 
totic behavior of observables at high energies. The expressions for the purely Coulomb 

background P matrix, its poles and residues, and purely Coulomb eigenfunctions of the 

> '. 

P-matrix approach are obtained. The nuclear-Coulomb low-energy scattering parame- 

O 

^ | the method of isolating the background P matrix. Simple explicit expressions for the 

nuclear-Coulomb scattering length and effective range in terms of the residual P matrix 
are derived. These expressions give a general form of the nuclear- Coulomb low-energy 
scattering parameters for models of finite-range strong interaction. Specific applications 



ters of two charged hadrons are investigated on the basis of this approach combined with 



of the general expressions derived in this study are exemplified by considering some ex- 



^ . actly solvable models of strong interaction containing hard core repulsion, and, for these 

H : 

models, the nuclear-Coulomb low-energy scattering parameters for arbitrary values of 
the orbital angular momentum are found explicitly. In particular, the nuclear-Coulomb 
scattering length and effective range are obtained explicitly for the boundary-condition 
model, the model of a hard-core delta-shell potential, the Margenau model, and the 
model of hard-core square-well potential. 



Electronic address: pet@gluk.org 



1 



1. INTRODUCTION 



The P-matrix approach to describing hadron-hadron interaction was first proposed 
by Jaffe and Low [1] and was then developed in a number of studies [2-4]. This approach 
is a modification of the well-known Wigner-Eisenbud 3?-matrix theory [5, 6]. Within the 
P-matrix approach, the scattering amplitude is expressed in terms of the logarithmic 
derivative of the wave function at the surface of the strong-interaction region — that 
is, in terms of the P matrix. Here it is assumed that the configuration space of the 
system is broken down into two regions, the external region, where the interaction can 
be described in terms of a two-particle potential, and the internal region, where strong 
interaction is dominant. For the P matrix, a so-called dispersion formula that appears 
to be its pole expansion and which establishes its energy dependence can be derived on 
the basis of quite general assumptions. Observables can then be described in terms of a 
finite number of parameters. 

In [3, 4, 7], a method was proposed for explicitly isolating a free background part in 
the P matrix. This method is advantageous in that it simplifies the implementation of 
the P-matrix approach in specific applications and extends the region of its applicability. 
The free P matrix, which corresponds to the absence of interaction, was isolated as 
the background part in the aforementioned studies. This is natural in dealing with the 
scattering of neutral particles. Here, we propose a generalization of the isolation method 
to the case involving charged particles, so that there is long-range Coulomb interaction 
in the system along with strong interaction. It is well known that, in this case, scattering 
theory requires a nontrivial modification. We show that the idea of explicitly isolating a 
background part in the P matrix as put forth in [3, 4, 7] can be implemented for charged 
particles as well and that, for the background P matrix, it is advisable in this case to 
take the purely Coulomb P matrix — that is, the logarithmic derivative of the regular 
Coulomb wave function at the surface of the interaction region. It turns out that the 
isolation of the background Coulomb part offers the same advantages as in the absence 
of Coulomb interaction. 

As an application of the P-matrix approach combined with the method for isolat- 
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ing the background Coulomb P matrix, we study the scattering length and the effective 
range for low-energy nuclear-Coulomb scattering. These parameters are important phys- 
ical quantities characterizing the scattering of charged hadrons and light nuclei at low 
energies. We obtain simple explicit expressions for the nuclear- Coulomb low-energy scat- 
tering parameters in terms of the parameters of the residual P matrix; these expressions 
make it possible to analyze and evaluate the nuclear- Coulomb scattering length and 
effective range and to find them directly for finite-range strong-interaction potentials. 
As a matter of fact, the expressions that we obtain here determine a general form of the 
nuclear-Coulomb low-energy scattering parameters for models of strong interaction of 
finite range. In [8], expressions for the nuclear- Coulomb low-energy scattering parame- 
ters in terms of the P-matrix parameters were obtained without resort to the isolation 
method. Those expressions are more cumbersome and less convenient in applications 
than the present ones. 

Much attention has been given to the nuclear-Coulomb low-energy scattering param- 
eters (see, for example, [9-12]) since these physical quantities play an important role in 
theoretical and experimental investigations. In some studies (see [13-17]), these param- 
eters were determined explicitly for some specific cases of separable nuclear potentials 
(in particular, for the Yamaguchi potential). Here, we find a general form of the nuclear- 
Coulomb scattering parameters for a rather broad class of local strong-interaction models 
- namely, for models of finite range. General expressions obtained for the low-energy 
parameters make it possible to determine these quantities explicitly for a number of ex- 
actly solvable strong-interaction models containing hard core repulsion. We emphasize 
that it is the use of the simple expressions obtained by the isolation method that made 
it possible to simplify significantly the relevant consideration. It should also be noted 
that the investigation presented here was performed for an arbitrary value of the orbital 
angular momentum i. 



3 



2. DISPERSION FORMULA FOR THE P MATRIX AND 
NUCLEAR-COULOMB OBSERVABLES 



For the elastic scattering of two charged strongly interacting particles, the radial 
wave function ipik (r) of relative motion in a state characterized by a specific value of the 
orbital angular momentum t is regular at the origin and satisfies the radial Schrodinger 
equation 

\ j 2 ia±x\ 

i>ik (r) = (1) 



^ + E Jl±H. v(r) 



dr 2 



with the potential 



V(r) = V s (r) + V c (r) , (2) 

where V s (r) is a short-range strong potential (by assumption, it has a finite range R) 
and V c (r) = 2£/c/r is an ordinary Coulomb potential. Here, £ is the Coulomb parameter 



= /ieie 2 



where a# is the Bohr radius, 



h 2 k ask ' 

h 2 



a B 



//eie 2 



(3) 



(4) 



ei and e 2 being the charges of the particles involved in the scattering process. We shall 
use the system of units where the reduced Planck constant and the doubled reduced 
mass are both equal to unity ( h — 2/j, — 1), so that the energy of the relative motion, 
E, is expressed in terms of the wave number k as E = k 2 . 

At infinity, the radial wave function satisfies the scattering boundary condition 



^Pik (r) 



[cos n (k) Fi (f , kr) + sin v { (k) G t (£, kr)] , 



(5) 



where ip lk (r) is the asymptotic wave function for the continuous spectrum, while 
Fi (£, kr) and Gi (£, kr) are, respectively, the regular and the irregular Coulomb wave 
functions [18], whose asymptotic behavior at infinity is given by 



In 

F t (£, kr) — > sin [kr — £ ln2A;r — — + o\ (k) ) , 



In 

Gi (£, kr) — > cos ( kr — £ In 2kr — — + a t (k) 



(6) 
(7) 
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Here, a t (k) = argT + l + is the purely Coulomb phase shift, and the total phase 
shift Si(k) has the form 

8 t (k) = a x (k) + v x (k) , (8) 

where v x (k) is the nuclear-Coulomb phase shift. 

The P matrix Pi(E) is defined in terms of the logarithmic derivative of the radial 
wave function at the surface of the strong-interaction region (r — R), 

P,(E) , M^f . (9) 

In the internal region r ^ R, we introduce a complete set of orthonormalized eigenfunc- 
tions U\ n (r) that satisfy the Schrodinger equation (1) and the homogeneous boundary 
conditions 

Uj„(0) = 0, u ln (R) = (10) 

at the ends of the interval [0, R\. Nontrivial solutions that obey the conditions in (10) 
exist only at some energy eigenvalues E\ n that are determined by solving the Sturm- 
Liouville problem specified by Eqs. (1) and (10). The orthonormalization conditions 
have the form 

R 

ui m ui n dr = 5 mn . (11) 


By expanding the wave function in the internal region in a series in eigenfunctions 
ui n (r), we find that, for the P matrix, there is the dispersion formula [1,3,4] 

00 771 2 

P^) = W) + ]T— — (12) 

n=l 

where 

lln = y/Ru' ln (R) . (13) 

Relation (12) involves a constant Pi(0), the P matrix at zero energy, and it is the iso- 
lation of this constant that ensures convergence of the remaining series. The dispersion 
formula (12), which represents a pole expansion of the P matrix, establishes the general 
form of its energy dependence. This dependence is completely determined by the states 
of the compound system which are characterized by the energy eigenvalues E\ n and the 



/ 
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residues jf n . These quantities in turn are controlled by the physical properties of the 
system in the internal region and are independent of energy E. 

Let us now establish the relation between the S matrix and the P matrix. For this, 
we note that the wave function in the external region, ip lk (r), can be represented in the 
general form 

^ifc (0 = \ [Ht ] & kr) - S l (k) Hj +) (£, kr)] ,r>R, (14) 

where (£, kr) are the Coulomb Jost solutions given by 

(£, kr) = [G t (£, kr) ± iF t (£, kr)} , (15) 

which represent the diverging and the converging waves distorted by the Coulomb po- 
tential. Accordingly, their asymptotic behavior is given by 

if (±) (£, kr) — e ±<Kr-Sln2kr-!?) 

By using the matching conditions at the point r = R and definition (9), we find that 
the S matrix can be expressed in terms of the P matrix as 

?(fc) m Pt } (k) - Pi (k) 
P/ +) (k) - P, (k) 

where 



Sl (k) = 3^ (k) ifc ^ ij ^ , (17) 



?(fc) , M = # <-> (e, fcg) 

is the S matrix corresponding to the scattering on a hard core of radius R in the 
presence of the Coulomb potential and P^ (k) are the logarithmic derivatives of the 
diverging and converging Coulomb waves at the boundary surface, 

l{±)nA _ kRH^' (£,kR) 
Hl ±] (£, kR) 

Hereafter, a prime denotes differentiation with respect to the variable p = kR. The real 

and the imaginary part of the function p/ + -* (k) are usually denoted by A; (k) and si (k), 

(+) * 

respectively; since the relation P{ (k) =P l (k) obviously holds, we can write 

P, (±) (*) = A, (*) ± is, (k) , (20) 
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where the functions si (k) and A/ (k) are expressed in terms of the Coulomb functions 
as [19] 



si (k) = 



kR 



if (e, kR) + of (e, ' 

A, (&) = (A;) |>, (£, fcfl) F; (£, kR) + G, (£, fcfl) G[ (£, kR) 



(21) 
(22) 



With the aid of Eqs. (17) and (20), it can easily be found that the nuclear-Coulomb 
phase shift can be represented as 

-W = c 'W + arctan P,w'-A,w (23) 

where the phase shift Q (k) for scattering on a hard core of radius R in the presence of 
the Coulomb interaction is given by 

^> s - arctan ftlt- < 24 > 

Expressions (17) and (23) for the observables reveal a significant drawback of the 
P-matrix approach based on the dispersion formula (12) as an approximation of the P 
matrix: if only a finite number of terms are retained, the observables in question will 
have an incorrect asymptotic behavior at high energies. By way of example, we indicate 
that, with increasing energy, the phase shift (23) will then behave as the phase shift 
(i (k) for scattering on a hard core; that is, it will tend to infinity, 

ui(k) — ► -kR + 0(l) . (25) 

But in fact, the phase shift must vanish at high energies, at least for regular potentials. 
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3. PURELY COULOMB P MATRIX 



An incorrect asymptotic behavior of observables at high energies can be avoided by 
isolating the background part in the P matrix. In the presence of Coulomb interaction, 
we define the background P matrix as the purely Coulomb P matrix — that is, in terms 
of the logarithmic derivative of the regular Coulomb wave function, 

We recall that the regular Coulomb wave function Fi (£, p) is expressed in terms of the 
confluent hypergeometric function $ (a, b; z) as [18] 

Fi (£, p) = Q (0 e- V +1 * (J + 1 - i& 2J + 2; 2tp) , (27) 

where 

C ^ = j^y^\r(l + l + i0\ (28) 

is the Coulomb penetrability factor that is a rather complicated function of energy and 
which is introduced in order to ensure the required asymptotic behavior of the function 
Fi(£,p) at infinity [see Eq. (6)]. If, however, we consider only the internal region 
< r < R, it is more convenient to introduce a solution that does not involve the factor 
C/(£) and which possesses simpler properties near the origin. We define such a solution 
4>i (£> P) through the relation 

F^,p) = (2/ + 1)HQ(0<M£,P) , (29) 

where the factor (2/ + 1)!! was introduced in order that, upon switching the Coulomb 
interaction off, the function <pi (£, p) reduce to the spherical Riccati-Bessel function: 

h (0, P) = 3i (P) ■ (30) 

The expression for the function 0; (£, p) in terms of a confluent hypergeometric func- 
tion can easily be found with the aid of Eqs. (27) and (29). The Coulomb P matrix as 
expressed in terms of the solution (pi (£, p) has the form 

(c) _ kR<j>\ (tj,kR) 

Pl {E) - MZ,kR) ' (31) 
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which is analogous to (26). 

The positions = kj^ of the poles of the Coulomb P matrix — they depend on 
the Bohr radius as and on the interaction range R (E^ = E^ (a,B,R)) — are defined 
by the roots of the denominator of the expression on the right-hand side of (31), 

<t>i(-^-MnR] =0. (32) 



The Coulomb eigenfunctions u?^ (r), which obey the Schrodinger equation (1) with 
the purely Coulomb potential and the boundary conditions (10), are given by 

J c) ( r \ - 7 '" <t>i (6n, k ln r) 

where ^ n = I/obA^. And the parameters 7^, which are determined from the normal- 
ization condition for the eigenfunctions, can be found if we use the Green's theorem, 

R 

u ik, (R) u'iki ( R ) - M k (R) u ik 2 ( R ) = { k i ~ k l) J u lkl u lk2 dr (34) 



for two solutions to Eq. (1) that correspond to two different energy values, k\ and k\. 
By substituting (33) into (34), going over to the limit k\ — > k 2 = k\ n and taking into 

(c) 2 

account Eqs. (11) and (13), we find that the expression for the Coulomb residues 7 in 
can be recast into the form 

7 M 2 = 2 Jl£_ (35) 

ln l-e l {i ln ,k ln R)/a B RE^<t> l {i ln MnR) 

Here the function 61 (£, p) is the derivative of the function <pi (£, p) with respect to the 

parameter £, 

ate \- d< ^ 1 (ia\ 
e i{£,P) = q£ , ( 36 ) 

and can be directly expressed in terms of a confluent hyper geometric function. 

Thus, we have completely determined the parameters of the Coulomb P matrix (its 
poles and residues) and found the Coulomb eigenfunctions. The dispersion formula for 
the Coulomb P matrix has the form (12); that is, 

p} c \E) = P, (c) (0) + V J'" , (37) 
' 1 E c) E — Ei n y ! 

n=l In 
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where the P matrix at zero energy, P^(0), is given by expression 

Jf(») = 7^-I, (38) 

which can be directly obtained from (31) for k — > by using the expansion of the regular 
Coulomb wave function and its derivative in terms of Bessel functions [18, 20, 21]. In 
expression (38), I v (z) are modified Bessel functions and the dimensionless parameter z 
is given by 

z = 2^2R/a B . (39) 

So far, we have considered the case of Coulomb repulsion (a B > 0). In the case of 
Coulomb attraction (a B < 0), P^°\o) has the form 

^'(o) = i^K)-;, (40) 



where J V {Q are Bessel functions and ( = 2^2R/ \a B \ . 

4. ISOLATING THE PURELY COULOMB BACKGROUND P MATRIX 

In nuclear-Coulomb P matrix (9), we now isolate explicitly the purely Coulomb 
background P matrix (26), following a way that is similar to that used to isolate explic- 
itly the free background P matrix in the absence of Coulomb interaction [3, 4, 7]. We 
represent this transformation in the form 

P(P) = P/ C) (P) + P(£). (41) 

With the aid of Eqs. (12) and (37), it can be shown that, for the residual nuclear- 
Coulomb P matrix P t (E), we have the expansion 



P l (E) = P l (0) + J2 



E ll E 7 <f 



E ln E-E ln E%E-E% 



(42) 



n=l 

By comparing the expansions (12) and (42) for the functions Pi(E) and Pi(E) we 
conclude that the isolation of the purely Coulomb background part in the nuclear- 
Coulomb P matrix according to (41) amounts to a partial summation of the series in 
(12), where one isolates the part that corresponds to the Coulomb interaction and which 
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is known explicitly. This naturally improves convergence of the original series, thereby 
making it possible to obtain a more accurate representation of observables as functions 
of energy. It can be shown that the expansion in (42) converges at the same rate as a 
series whose general term is proportional to 1/n 4 , while the expansion in (12) converges 
as 1/n 2 — that is, much more slowly. 

By making transformation (41) in Eq. (17) and using the relation 



p/ c) (k) - (k) = k -^-^ = -m— > ( 43 ) 

F l (£, kR) (£, kR) <f>(fc) 



which can easily be verified for logarithmic derivatives, we can straightforwardly express 
the nuclear-Coulomb S matrix (E) = e 2lu '^ in terms of the residual P matrix 
Pi{E) as 



:(cs) ( ^ _ i+ c ;- } (k)Pj(E) 

1 + 4 +) (k) P^E) 

where the functions c| ± ' ) (k) are determined according to (43) and obviously satisfy the 
relation c\ (k) =c { (k). With the aid of (44), we can easily represent the nuclear- 
Coulomb phase shift as 



F?(Z,kR) Pi(E) 

Vl (k) = - arctan ^ 1 y ' — ^ . 45 

kR + F l ^,kR)G l ^,kR) P t {E) 

From Eqs. (44) and (45), it is obvious that, if only a finite number of terms are 
retained in the pole expansion (42) for Pi(E), the S matrix and the phase shift will 
have a correct asymptotic behavior at high energies: 

S t (k) — > 1 , v x (k) — > . (46) 

Thus, an isolation of the purely Coulomb background term in the P matrix leads to a 
correct asymptotic behavior of the observables at high energies if the residual P matrix 
is approximated by a finite number of pole terms. We can see that the transformation 
(41) has a transparent mathematical and physical substantiation and that its application 
provides the same advantages as in the absence of Coulomb interaction. 
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5. EXPRESSIONS FOR THE NUCLEAR-COULOMB LOW-ENERGY 
SCATTERING PARAMETERS IN TERMS OF THE RESIDUAL P 

MATRIX PARAMETERS 



A great number of studies have been devoted to the problem of generalizing and mod- 
ifying effective-range theory in the presence of long-range Coulomb interaction (see, for 
example, [22-27]). As a result, the Coulomb-modified effective-range function K csl (E) 
was introduced, and the nuclear- Coulomb scattering length and effective range were 
determined for the case of S'-wave scattering, as well as for scattering in a state char- 
acterized by an arbitrary value of the orbital angular momentum £. In [22, 25], it was 
shown that, in the case of an arbitrary orbital angular momentum, the effective-range 
expansion in the presence of Coulomb interaction has the form 



Ki(E) = (2Z + 1)!! 2 C?(0 k 2l+1 

1 1 



2£ 

cot ^ + cf(0 h ® 



= -- + l ri k 2 + ... : (47) 
0,1 * 

where the function h(£) is expressed in terms of the digamma function ip(z) = r'(z) /T(z) 
as 

h(Z) = Rei/>(l + i£)-hi\Z\. (48) 

In the complex plane of energy E, the nuclear- Coulomb effective-range function 
Ki(E) is analytic in some domain near the origin [26]; hence, it can be expanded in 
the Maclaurin series (47) in powers of E in the vicinity of the point E — 0. Thus, a 
special role of the function Ki(E) is associated with its analyticity near E — 0. The 
nuclear-Coulomb scattering length a\ and effective range r t are determined in this case 
in terms of the coefficients in the expansion (47) of the function Ki(E) . We note 
that, in a large number of studies, the nuclear-Coulomb quantities, such as ai, ri, Ki 
and others are equipped with the additional indices c and s, which label parameters 
associated with, respectively, Coulomb and short-range interactions. This was done in 
order to distinguish these quantities from their counterparts in the absence of Coulomb 
interaction, which are labeled only with the index s. Since we do not consider here 
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the case where there is no Coulomb field, the indices c and s are suppressed on all 
nuclear-Coulomb quantities. 

For a further analysis, it is reasonable to introduce the dimensionless inverse scat- 
tering length 7/, and the dimensionless effective range pi as 

= (49) 
Pi^ZlPa^n. (50) 

In the particular case of scattering by a hard core of radius R (ipi k (R) — , Pi (E) = oo), 
the nuclear- Coulomb low-energy scattering parameters 72 and pi can easily be found in 
the explicit form [8] 

(51) 

rf = l-^ + ^, (52) 



V 



where the superscript h denotes, as previously, a hard core and I v (z) and K u (z) are 
modified Bessel functions. The constants v, \i, and p\ are given by 

v = 21 + 1 , (53) 

A, = 1(1 + 1) , (54) 
m = 4u\i , (55) 
while the dimensionless parameter (5 is defined as 

P = —. (56) 
a B 

As before, the parameter z has the form (39). In order to render the expressions pre- 
sented below less cumbersome, it is convenient to isolate explicitly, in the low-energy 
parameters 7? and pi, the parts that correspond to scattering on a hard core. Accord- 
ingly, we set 

7/ = li + li , (57) 
Pi = Pi + Pi, (58) 
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defining, in this way, the residual low-energy nuclear- Coulomb scattering parameters 7; 
and pi. 

Let us further express the nuclear-Coulomb scattering length and effective range in 
terms of the residual P matrix. By substituting (45) into (47), we find that the nuclear- 
Coulomb effective-range function as expressed in terms of the residual P matrix is given 
by 



K l (E) = (2l + l)\\ 2 Cf(0k 



(59) 



where, as before, we use the notation p = kR. Let us expand the right-hand side of Eq. 
(59) in a Maclaurin series in powers of energy E = k 2 . It is obvious that, as long as 
we are interested neither in the shape parameter nor in higher expansion coefficients, it 
is sufficient to retain only the terms that are linear in E. We further make use of the 
known relation for the Coulomb penetration factor [18], 

c'KO 2 2 ' 



^)-(2, + l)!^ + « 2 »( 22 + «-( ii + «' 
and of the asymptotic expression for the function h(£) at low energies [18], 



(60) 



+ 



+ 



(61) 



12£ 2 120£ 4 

The expansions of the Coulomb wave functions in power series in energy E were previ- 
ously studied by many authors [20, 21, 28, 29]. To terms that are linear in energy, these 
expansions for the case of Coulomb repulsion can be written as 



(2l + l)\Q(Q z , 

( 2 £)i+i 2 S A ) 



z 



(2/ + l)!Q(Q 
(20'Cg(0 



z<K v (z) + 



see 2 

,3 



Z 



96f 2 



21 

I v +i{z)-\ I v + 2 {z) +... 

z 

21 

K v+1 (z) K v+2 (z) +. 



(62) 



(63) 



The dispersion relation for the residual P matrix (42) can be recast into the form 



n=l 



(c) 2 



E-E, 



In 



E-E, 



.(c) 
In 



(64) 



This expansion contains no additional parameters and is completely determined by the 
quantities Ei n and jf n . It can be shown that the expansion in (64) converges at the 
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same rate as a series whose general term is proportional to 1/n 2 . The analyticity of the 
residual P matrix in a vicinity of the point E = immediately follows from (64) if all 
energy eigenvalues differ from zero. The expansion of the residual P matrix in a power 
series in energy E can be written in the form 



P l (E) = P l + Q lP 2 + ... , 



(65) 



where 



Pi = Pi(0), 
1 



(66) 

Qi = ^Pl(0) (67) 
are dimensionless expansion coefficients. We can easily express the quantities Pi and 



Qi in terms of the P-matrix parameters 

oo 



as 



7ta 



E 



M E 



n=l 



'In 

(c? 
Tin 



E. 



'In 



El 



(68) 



(69) 



71=1 L^to -I 

We also note that, on the basis of Eq. (41), the quantities Pi and Qi can be determined 
from the relations 

(70) 
(71) 

(72) 
(73) 

(74) 



where 



Pi = P/ c) + Pi , 

Qi = Q\ c) + Qi , 

Pi = Pi(0), 

l 



are parameters in the expansion of the P matrix, 

Pi(E) = Pi + Qip 2 + ... . 

The parameters P^ and in the expansion of the purely Coulomb background P 
matrix can be found explicitly from Eqs. (26) and (62). The results 



are 



p (c) _ , 1 z (z) 



(75) 
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2(3 (3Q^ + l)=lu+ (iz - ^±lM + 2 09 - Aj) . (76) 

Substituting now expressions (60)- (63) and (65) into Eq. (59) and taking into ac- 
count Eq. (47), we arrive at explicit expressions for the inverse scattering length and 
the effective range pi, also referred to as the dimensionless nuclear- Coulomb residual 
parameters. The results are given by 

li = , (77) 

k =3/3S + 4^^i-/, (78) 



z 2 ii P l z I v 

where I u = I v (z). In the particular case of interaction in the S state ( / = 0), the last 

formulas are somewhat simplified to become 

7 = — ^ , (79) 
21* P 

% = W%-~- (80) 



z 2 7 r p 2h 

In the case of Coulomb attraction (as < 0), all the above formulas are valid upon the 
substitution of conventional Bessel functions for modified ones: 

h{z) - ij„(C), K v (z) - ^n(C), C = 2vWN. (81) 

Formulas (77) and (78) yield general explicit expressions for the nuclear-Coulomb 
low-energy scattering parameters in terms of the residual P matrix parameters. These 
expressions make it possible to obtain directly a general form of the nuclear-Coulomb 
scattering length and effective range for models of finite-range strong interaction. 

6. NUCLEAR-COULOMB LOW-ENERGY SCATTERING 
PARAMETERS FOR EXACTLY SOLVABLE MODELS CONTAINING 

HARD CORE REPULSION 

For specific applications of the above general expressions, we will consider some 
exactly solvable models of finite-range strong interaction containing hard core repulsion. 
For these, we will find explicitly the nuclear-Coulomb low-energy scattering parameters 
for arbitrary values of the orbital angular momentum. 
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6.1. Boundary Condition Model 

In the boundary condition model, the interaction in the internal region is determined 
by a single energy-independent parameter, the value of the logarithmic derivative of the 
wave function at the boundary surface — that is, the constant P\. It is obvious that the 
parameter Qi vanishes in this case. Thus, we have 

Pi (E) = Pi , (82) 

Qi = 0, (83) 

whence it follows that 

P l = P l - P} c) , (84) 

Qi = -Q\ C) ■ (85) 
The nuclear-Coulomb low-energy scattering parameters can then be written as 

^ = T~~ w\ ' (86) 

2/, 2 (Pi-P^) 

4r = W , f . +4^1^-1, (87) 

where the quantities Pf 1 and Q^f 1 are given by Eqs. (75) and (76). 

6.2. Hard- Core Delta-Shell Potential 

Let us consider the case where the strong interaction is described by the delta-shell 
potential concentrated on the sphere of radius R and supplemented with a hard core of 
radius R c less than R, 

{+oo , r < R c 
' • (88) 

-±5(r-R) ,r>R c 

Here, A is the dimensionless interaction constant. In this case, the wave function in the 
internal region ( r < R) is a linear combination of the Coulomb wave functions, 

ij lk (r) = At (k) [Gi (£, kR c ) Ft (£, kr) - F t (^, kR c ) G t (^, kr)] , R c < r < R (89) 



17 



and satisfies the zero boundary condition at r = R c : ipi k (R c ) = . At the bound- 
ary surface (r = R), the wave function is continuous, but its derivative undergoes a 
discontinuity, 

i,[ k (R + 0) - V4 (i2 - 0) = ^, fc . (90) 

By using formulas (89) and (90), we find that, in the case of a hard-core delta-shell 
potential, the P matrix can be represented as 

^ P G l (Z,x)F l (Z,p)-F l (Z,x)G l (Z,p) ' [ } 

where x = kR c . With the aid of the definition of the background Coulomb P matrix 
(26) and the representation in (41), we obtain the residual P matrix for the potential 
(88) in the form 

p (F) - pFi{^x) /F,(£,p) 

l{ ' G l (Z,x)F l (Z,p)-F l (Z,x)G l (Z,p) • 1 ' 

By means of the expansion of the Coulomb functions (62) and (63), we derive the 

parameters of the residual P matrix, Pi and Qi, according to (65). The results are 

P, = - h {y) /L (Z) - A (93) 

2 K v (y) I u (z) — I u (y) K v (z) ' 1 ) 

2A(3 2 Q l = (lz 2 - fn) (\ + P,) + (A, -0)x 

(A, - a) I (A, -0)-2z \K V (y) I v+1 (z) + I u (y) K u+1 (z)] I v (y) [h (z) + I 2 U (y) /I 2 (z) 

\K V (y) I u (z) - I v (y) K v (z)} 2 

(94) 

where a = R c / cle and y = 2y / 2a. By substituting (93) and (94) into (77) and (78), 
we now find that the nuclear-Coulomb scattering length and effective range for the 
hard-core delta-shell potential can be represented as 

^ = ts^m (95) 

% = r il V) ~il Z L + 2A*/„ (z) [lzl v (z) + 4 (P - A,) I v+1 (z)} , (96) 

where 7 ; (y) and pf (y) are the low-energy nuclear-Coulomb parameters for a hard core 
of radius R c (Eqs. (51) and (52) with substitutions z — > y and f3 — > a), while 7/ 1 (z) = 7/ 1 
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and pi (z) = Pi are the parameters for the hard core of radius R (Eqs. (51), (52)). In the 
limiting case of a delta-shell potential without a core (R c — > 0, y — > 0, and 7^ (y) — > oo), 
expressions (95) and (96) reduce to the known expressions for the low-energy nuclear- 
Coulomb parameters for scattering on a delta potential [16, 17] 

% = -2XJUz) ' (97) 
|r = 2XzI v (z) [lzl v (z) +4 ((3- A/) I v+1 (z)} . (98) 

6.3. Margenau Model 

In the Margenau model [30], strong interaction is simulated by a square- well potential 
containing hard core repulsion; in addition Coulomb interaction is assumed to be absent 
in the internal region. The latter is justified by the fact that, in the internal region, the 
Coulomb interaction is much weaker than strong interaction. Thus, the total interaction 
in this model is described by the potential 



V(r) 



+oo , r < R c , 

-V ,R c <r<R, (99) 
2£k/r , r>R. 



In this case, the wave function in the internal region has the form 

1> lk (r) = Ax (k) [m (KR C ) ji (Kr) - 3l (KR C ) m (Kr)} , R c < r < R , (100) 

where K = y/Vo + E, and ji and ni are the spherical Riccati-Bessel functions. The P 
matrix can then be written as 

M = P ^|^, (ioi) 

Ji {x) n t (p) - n { (x) 31 (p) 

where x = KR C and p = KR. From the above, we can easily determine the P-matrix 
parameters Pi and Q\. The results are 

p ji (xojn'j (p )-ni (x )j ; / (po) 
Ji {x )ni (p ) - ni {x )ji (p ) 
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2Q l + l = X ' + P '- P ' + " (103) 

Po [Ji {x )ni (p ) - m {xo) j i (po)J 

where 6 = R c /R, x = K R C , p = KqR, and i^o = \fVo- By using Eqs. (70) and 

(71) and taking into account Eqs. (75) and (76), we can find the parameters of the 

residual P matrix — P\ and Qi; the nuclear-Coulomb low-energy scattering parameters 

in this case are given by 

1 0n Ji (xp) n l+1 (po) - m (x ) jj+i (po) I v+1 

= ~ 2 Po — —, — s 7 — x 7 — r— ; — v z—— , (104) 



11% J h (x ) ni (po) - n t (x ) j t (p ) I v 
Pi _ W +2 {(3 _ Xi) S^ +3(3 b_ - [ 3l (x ) n l+l (po) - m (x ) j l+l (p )] 2 + 



z2l lll z I* " II \ji{xo)ni(p ) -ni(x )ji(p )f 



+ 



3z/£ + 2/po + (P - \i) 



Z Iu+1 



ji (x ) n l+1 (po) - ni (x ) ji+i (po) 
ji (x ) ni (p ) - ni (x ) ji (p ) 



£ Iv 
Where ^=^K - 

6.4- Hard-Core Square- Well Potential 

For the case of a strong-interaction potential in the form of a square well with a hard 
core, 

{+oo , r < R c , 
(106) 
-V 9(R-r) ,r>R c 

we confine ourselves to determining the nuclear-Coulomb scattering length. For the 

simpler case of a square-well potential without a core, the nuclear-Coulomb scattering 

length and the nuclear- Coulomb effective range were found in [8]. In this case, the wave 

function in the internal region has the form 

Vfc (r) = A l (k) [G t ( ~ , KR C ) F l (S, Kr) - F t (S, XP C ) ( ~ , Kr)] , i? c < r < i* , 

(107) 

where S = , and the P matrix is given by 

„, F1 fj (5, x) gj (5, g) - G, (5. g) g (5, p) 

where, as in the preceding subsection, x = KR C and p = KR. In accordance with (70) 

and (75), the parameter Pi of the residual P matrix then assumes the form 

p = Fd^x )G' l ^ ,po)-Gi^ ,x )F;^ ,p ) z l v+1 (z) 
1 Po Fi (£ , x ) Gi (Co, Po) - Gi (Co, x ) (Co, a,) 2 J„ (z) 1 + J ' 1 J 
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while the nuclear-Coulomb scattering length is given by 



1 



2p 



F l (fr, gp) G; (Co, go) - g; (6, go) (Co, go) 
(6, &o) (f , Po) - G/ (£ , xo) ^ (fo, Po) 



— z 




+ (HO) 



7. CONCLUSION 



In summary, an explicit isolation of the purely Coulomb background part in the P 
matrix leads to a correct asymptotic behavior of physical observables at high energies 
when the residual P matrix is approximated by a finite number of pole terms. Concur- 
rently, the isolation of the background P matrix makes it possible to improve convergence 
of the remaining expansions. The transformation in (41) has a transparent mathematical 
and physical substantiation, and its application provides the same advantages as in the 
absence of Coulomb interaction. In addition, the explicit isolation of the purely Coulomb 
background part in the P matrix makes it possible to obtain the simple general expres- 
sions (77) and (78) for the nuclear- Coulomb low-energy scattering parameters in terms 
of the residual P matrix. With the aid of these expressions, we can directly calculate the 
nuclear-Coulomb scattering length and effective range for finite-range strong-interaction 
potentials. If the Schrodinger equation for these potentials admits of an exact solution 
in the presence of Coulomb interaction, the nuclear-Coulomb parameters can be found 
explicitly. In general, the nuclear- Coulomb low-energy scattering parameters can be ob- 
tained for arbitrary short-range strong-interaction potentials at any value of the orbital 
angular momentum i. On the basis of the expressions derived in the present study, we 
have found explicitly the nuclear-Coulomb scattering length and effective range for the 
boundary-condition model, for the model of a hard-core delta-shell potential, for the 
Margenau model, and for the hard-core square-well potential at arbitrary values of the 
orbital angular momentum. 
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